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Abstract 

During September 2010, an outbreak of acute hemorrhagic conjunctivitis reemerged in Jiangsu, three years after the 
nationwide epidemic in China in 2007. In total, 2409 cases were reported, 21 18 of which were reported in September; 79.8% 
of those affected were students or teachers, with a median age of 16 years. To identify and demonstrate the genetic 
characteristics of the etiological agent, 52 conjunctival swabs were randomly collected from four different cities. After 
detection and isolation, 43 patients were positive for coxsackievirus A24 variant according to PCR and 20 according to 
culture isolation. Neither adenovirus nor EV70 was detected. A phylogenetic study of the complete 3Cpro and VP1 regions 
showed that the Jiangsu isolates clustered into a new lineage, GIV-C5, with two uniform amino-acid mutations that 
distinguished them from all previous strains. Another new cluster, GIV-C4, formed by Indian isolates from 2007 and Brazilian 
isolates from 2009, was also identified in this study. Interestingly, our isolates shared greatest homology with the GIV-C4 
strains, not with the isolates that were responsible for the nationwide acute hemorrhagic conjunctivitis epidemic in China in 
2007. Although all our isolates were closely related, they could be differentiated into two subclusters within GIV-C5. In 
conclusion, our study suggests that a new cluster of coxsackievirus A24 variant that had already evolved into diverse strains 
was associated with the acute hemorrhagic conjunctivitis outbreaks in Jiangsu in September 201 0. These viruses might have 
originated from the virus isolated in India in 2007, rather than from the epidemic strains isolated in China in 2007. 
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Introduction 

Acute hemorrhagic conjunctivitis (AHC), characterized by the 
sudden onset of ocular pain, swelling of tlie eyeUds, a "foreign 
body" sensation or irritation, epiphora (excessive tearing), eye 
discharge, and photophobia, is mainly caused by enterovirus 70 
(EV70), coxsackievirus A24 variant (CA24v), or adenoviruses [1- 
3]. In China, the first AHC outbreak was reported in 1 97 1 [4]. 
Because AHC is a notifiable infectious disease in China, all cases 
diagnosed by physicians have been registered in the National 
Disease Supervision Information Management System (NDSIMS). 
In 2007, a nationwide AHC epidemic caused by CA24v was 
reported in China, with a total of 74,263 AHC cases [5,6], after 
which the number of AHC cases returned to baseline in 2008 and 
2009. In 2010, outbreaks of AHC reemerged around the world [7- 
11]. In China, outbreaks were reported in Zhejiang, Guangdong, 
Guangxi, Shandong, Henan, Fujian, Beijing, and Chongqing 
[7,12]. Several studies demonstrated that the outbreaks were 
caused by CA24v. However, to the best of our knowledge, in 
China in2010, the complete sequences of the VPl and 3Cpro 
regions of CA24v, which are usually used for phylogenetic 
analyses, were only determined in Zhejiang and Guangdong 
Provinces and submitted to GenBank. These sequences indicated 
that the CA24v strains causing these outbreaks belonged to Group 
IV [7,8]. The objectives of the present study were to investigate the 



AHC outbreak in Jiangsu Province, China, in September 20 1 0 
and to identify the etiological agent causing this outbreak and 
determine its genetic characteristics. 

Results 

The Outbreak 

Of the 4276 cases of AHC in Jiangsu Province reported to the 
NDSIMS between January 1 2007 and December 31 2010, 2409 
were officially reported in 2010, which was approximately 25-fold 
higher than in 2009 and 20-fold higher than in 2008, and even 
0.5-fold higher than in 2007 (Table 1). Unlike previous years, 
approximately 88% of the patients in 2010 were reported during 
September in 13 cities of Jiangsu Province. Of all the patients, 
79.8% were students, kindergarten children, or teachers (Table 1). 
Reports began to increase from September 5, reached a peak on 
September 1 3, and returned to baseline on October 1 (Fig. 1). All 
the patients had conjunctival congestion and a clear history of 
contact with people presenting similar symptoms. No death from 
AHC was reported. 

Virus Detection and Isolation 

Neither adenovirus nor EV70 was detected in 52 patients with 
acute conjunctivitis. 
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Figure 1. Distribution of tKie reported AHC cases in Jiangsu Province during the outbrealc of September 2010. 

doi:1 0.1 371 /journal.pone.0086883.g001 



However, 82.7% (43 of 52) of the patients were positive for 
CA24v according to PGR and 38.5% (20 of 52) according to 
culture isolation, suggesting that CA24v was the causative agent of 
the outbreak. 

Phylogenetic Analysis of the 3Cpro Region 

As in previous studies [8,13], phylogenetic analysis of the 3Cpro 
region showed that all the CA24v strains could be divided into 
four major groups. Our strains belonged to Group IV and formed 
a new cluster, G5, with strains isolated throughout the world after 
2010. Another new cluster (G4) formed by strains from India 
(2007) and Brazil (2009) was also identified in this study (Fig. 2). 



Analysis of the sequence identities of all five clusters of Group IV 
(Table 2) showed that our isolates shared the greatest homology 
with the C4 strains, rather than with the G3 strains, which were 
responsible for the nationwide AHC epidemic in China in 2007. 

Further analyses showed that our strains clustered in two 
subclusters: C5A and C5B (Fig. 2). C5A, consisting of isolates 
P19, P29, and P41, with sequence identities of 99.6%-100%, 
shared greatest homology with strain GD03/2010 (99.6%). C5B, 
formed by other strains sharing 99.3%-100% homology, clustered 
with strain 378/India/2010 (99.3%-99.5%). The homology 
between our strains and Zhejiang strains isolated in 2010 ranks 



Table 1. AHC cases reported to NDSIMS in Jiangsu from 2007 to 2010. 



Year 


Total cases 


Cases in September 


Student cases 


Median age 


MaleiFemale 


2007 


1649 


984 


702 


25 years 


1.99:1 


2008 


122 


31 


38 


30 years 


1.44:1 


2009 


96 


8 


20 


27 years 


2.10:1 


2010 


2409 


2118 


1690 


1 6 years 


1.69:1 
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Figure 2. Phylogenetic analysis based on tlie 549-nucleotide 3Cpro gene of CA24v. All Jiangsu CA24v isolates identified in this study, 
marl<ed with blacl< triangles, were compared with strains available in GenBank. Strains for which the complete sequences of both the 3Cpro and VP1 
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regions were available are indicated with **. The MEGA 4.0 software was used for the phylogenetic analysis. The stability of the nodes was assessed 
using neighbor-joining cluster analysis with 1000 bootstrap replications, and only bootstrap values >70% are shown at the nodes. 
doi:1 0.1 371 /journal.pone.0086883.g002 



second to the Indian strain (98.9%-99.1%), which is identical to 
the homology between C5A and C5B (98.9%-99.1%). 

Phylogenetic Analyses of the VP1 Region 

The sequence homology of the VPl region confirmed the 
findings of the 3Cpro region analysis (Fig. 3 and Table 3). The 

C5 lineage, which included our strains, also shared the highest 
sequence identities with the C4 lineage, which contained the same 
strains as in the 3C region analysis. The subclusters C5A and C5B 
were similarly identified with the VPl region. The sequence 
identity between C5A and C5B ranged from 97.4% to 98.3%. As 
before, C5A shared greatest homology with strain GD03/2010 
(99.2%-99.7%). However, because there were no complete VPl 
sequences for the Indian strains isolated in 2010, we demonstrated 
that C5B clustered and shared its greatest homology with strain 
Zhejiang/12/10 (99.0-99.6%). 

Amino Acid Sequence Analyses of the 3Cpro and VPl 
Regions 

An amino acid sequence analysis revealed two uniform 
mutations in all the strains that formed the C5 lineage after 
2010. One was the isoleucine at position 54 of the 3Cpro region, and 
the other was the aspartic acid at position 301 of the VPl region, 
whereas these two positions contained valine and asparagine, 
respectively, in all the other strains available in the GenBank 
before 2010 (Table 4 and Table 5). 

Discussion 

Our study demonstrates that CA24v was the causative agent of 
the AHC outbreak in Jiangsu Province in September 2010. CA24v 
is a major etiological agent of AHC. After its first description in 
Singapore in 1970, AHC outbreaks caused by CA24v have 
occurred periodically throughout the world [2,4,8-11,13-17]. 
According to previous studies, four genotypes of CA24v can be 
distinguished and identified by phylogenetic analysis of the 3Cpro 
and VPl regions [18,19]. Pei-Yu Chua reported that the CA24v 
strains of the fourth genotype, GIV, isolated from 2000 to 2007, 
can be divided into three different clusters: GIV-Cl, GIV-C2, and 
GIV-C3 [13]. The Chinese isolates isolated in Guangdong, 
Yunnan, and Zhejiang provinces in 2007 and 2008 were all 
clustered within GIV-C3 [5-6,8]. In this study, two new clusters 
were identified containing strains isolated after 2007. Our strains 
belonged to GIV-C5, and the strains from India (2007) and Brazil 

Table 2. Homology of the complete 3Cpro regions of the 



strains in different clusters within CA24v genotype GIV 
investigated in this study. 




Homology 

(%) 


Gives GIV-e4 


Gives Giv-e2 Giv-ei 


Gives 


98.0-100.0 - 




GIV-C4 


95.8-97.0 98.7 




GIV-C3 


92.8-94.2 95.8-96.8 


98.0-99.8 - 


GIV-C2 


93.8-96.8 95.4-97.0 


94.6-96.4 98.0-99.6 - 


GIV-Cl 


91.0-91.9 92.5-93.4 


91.9-93.6 93.0-94.2 99.5 
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(2009) formed GIV-C4 (Fig. 2 and Fig. 3). Contrary to previous 
studies [7,8], the differentiation of GIV-C5 and GIV-C4 is 
supported by phylogenetic analyses. First, our strains shared the 
highest sequence identities with the other strains within their 
clusters, whereas the identities were much lower with isolates that 
were categorized outside their clusters. The sequence identities 
were similar to the scores for the strains of other clusters (CI, C2, 
and C3) when compared with intra- and extra-clusters (Table 2 
and Table 3). Second, amino acid analyses revealed two uniform 
mutations that were only present in all the strains that formed the 
C5 lineage (Table 4 and Table 5). These results suggest that the 
CA24v strain that caused the AHC outbreak in Jiangsu in 
September 2010 belonged to the new C5 cluster. 

Interestingly, the new C5 cluster shared greatest homology and 
clustered with the C4 strains, which were represented by Indian 
strain MHO isolated in 2007 and Brazilian strain BR-PE/10 
isolated in 2009 [10,15], and not with C3, which was responsible 
for the nationwide AHC epidemic in China in 2007. However, 
other Indian strains, such as Ml 33 of GIV-C3, were also isolated 
in 2007. These results suggest that clues to the differentiation of C3 
and C4 can be found retroactively as early as 2007 in India. C5 
might then have evolved from C4, and has since evolved 
differently from the Chinese isolates of 2007 (Fig. 2 and Fig. 3). 

Our isolates were also differentiated as two subclusters within 
C5. The C5A cluster, formed by P19, P29, and P41, grouped 
together with the strains isolated in Guangdong Province in China 
in 2010 [8], whereas the others that formed C5B clustered with the 
strains isolated in Zhejiang Province in China in 2010. In an 
analysis of all the strains available in GenBank, we found that no 
strain isolated in 2010 in Guangdong clustered with C5B and no 
strain isolated in the same year in Zhejiang clustered with C5A. In 
other words, in 2010, at least when and where the studies were 
performed, the CA24v outbreaks in Guangdong might have been 
caused by strains that clustered with C5A alone, whereas in 
Zhejiang, CA24v only clustered with C5B [7,8]. However, in 
Jiangsu, the two subclusters were epidemic at the same time. 

C5A shared the greatest sequence identity in both the 3Cpro 
and VPl regions with strains from Guangdong Province in 2010 
[8]. The VPl region of C5B also shared the greatest identity with 
the strains isolated in Zhejiang Province in 2010 [7]. However, the 
3Cpro region of C5B shared its highest homology with the strain 
isolated in India in the same year. Although the complete VPl 
region of the Indian isolate was not available in GenBank, we 
deduced that frequent mutations and/or recombination events 
have occurred within the C5 cluster, generating its diversity. 

In conclusion, our study suggests that a new cluster of 
coxsackievirus A24 variant, which had already evolved into 
diverse strains, was associated with the acute hemorrhagic 
conjunctivitis outbreak in Jiangsu in September 2010. These 
viruses might have originated from the Indian strains isolated in 
2007, rather than from the epidemic strains present in China in 
2007. 

Materials and Methods 

Ethics Statement 

This study was approved by the Ethics Committee of the 
Jiangsu Provincial Center for Disease Control and Prevention. 
Written informed consent for the use of the 52 clinical specimens 
for virus isolation and detection was obtained from all the patients 
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Figure 3. Phylogenetic analysis based on the 91 5-nucleotide VP1 gene of CA24v. All Jiangsu CA24v isolates identified in this study, marked 
with black triangles, were compared with strains available in GenBank. Strains for which the complete sequences of both the 3Cpro and VP1 regions 
were available are indicated **. The MEGA 4.0 software was used for the phylogenetic analysis. The stability of the nodes was assessed using 
neighbor-joining cluster analysis with 1000 bootstrap replications, and only bootstrap values >70% are shown at the nodes. 
doi:1 0.1 371 /journal.pone.0086883.g003 



Table 3. Homology of the complete VPl regions of the 
strains in different clusters within CA24v genotype GIV in this 
study. 





Homology (%) 


Gives 


GIV-C4 


GIV-C3 


GIV-C2 


Gives 


97.3-99.9 








GIV-C4 


96.5-97.7 


98.9 






GIV-C3 


94.1-97.3 


96.1-97.5 


97.4-99.8 




GIV-C2 


94.1-95.4 


94.5-95.8 


94.5-96.3 


98.0-99.6 
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or the legal guardians of child patients. The Ethics Committee 
waived the requirement for the consent of the patients or 
guardians to use the basic information about the 4276 cases of 
AHC obtained from NDSIMS, because no samples were acquired 
from those patients. Patient identifiers, including names, addresses, 
and case numbers, were removed from any parts of the documents 
used in this study to ensure patient confidentiality. 

Clinical Information and Specimen Collection 

The 4276 cases of AHC in Jiangsu Province that were reported 
to NDSIMS between January 1 2007 and December 31 2010 were 
investigated in this study. In total, 52 conjunctival swabs from 
AHC patients aged 7-56 years (28 males and 24 females in the 
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Table 4. Amino acid mutations in the 3Cpro region in different genotypes of CA24v. 
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*Amino acid mutation only present in all GIV-C5 strains. 
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Table 5. Amino acid mutations in the VP! region in different genotypes of CA24v. 
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Y 


T 


T 




GIV-C2 


gi|58396775|Ningbo3-02 


T 


P 


L 


A 


V 


1 D 


R 


A 


H 




M 


Y 


T 


T 






gi|58396777|Hangzhoul3-02 


T 


P 


L 


A 


V 


1 D 


R 


A 


H 




M 


Y 


V T 


T 






gi|3263791 79|BR-RJ/37/2003 


T 


P 


L 


A 


V 


1 D 


R 


A 


H 




M 


Y 


T 


T 






gi|l 57885039|Spain41 92EV04 


T 


P 


L 


A 


V 


1 D 


R 


A 


H 




M 


Y 


T 


T 






gi|326379209|BR-ES/93/2005 


T 


P 


L 


A 


V 


1 D 


R 


A 


H 




M 


Y 


T 


T 




GIV-C3 


gi|347309220|GD332/2007 


T 


H 


L 


A 


V 


1 D 


R 


A 


H 




M 


Y 


T 


T 






gi|32000051 3|2007-M1 33-IND 


T 


H 


L 


A 


V 


1 D 


R 


A 


H 




M 


Y 


T 


T 






gi|258547237|Zhejlang/26/08 


T 


H 


L 


A 


V 


1 D 


R 


A 


H 


Q 


M 


Y 


T 


T 




GIV-C4 


gi|320000499|2007-Ml 10-IND 


T 


H 


L 


A 


V 


1 D 


R 


A 


H 




M 


Y 


T 


T 






gi|32637921 9|BR-PE/9/2009 


T 


H 


L 


A 


V 


1 D 


R 


A 


H 




M 


Y 


T 


T 




Gives 


gi|33877701 0|Zhejlang/l 2/1 0 


T 


H 


L 


A 


V 


1 D 


R 


A 


H 




M 


Y 


T 


T 


D* 




gi|347309252|CA24-GD03/201 0 


T 


H 


L 


A 


V 


1 D 


R 


A 


H 




M 


Y 


T 


T 


D* 




gi|375173647|P29/Jlangsu/10 


T 


H 


L 


A 


V 


1 D 


R 


A 


H 




M 


Y 


T 


T 


D* 




gi|471434775|lshlgakl/29/201 1 


T 


H 




A 


V 


1 D 


R 


A 


H 




M 


Y 


T 


T 


D* 



*Amlno acid mutation only present In all GIV-C5 strains. 
dol:10.1371/journal.pone.0086883.t005 
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acute phase), were randomly collected from four hospitals in 
different cities (13 from Nanjing, eight from Taizhou, 15 from 
Yancheng, and 16 from Lianyungang) during September 2010 (1— 

5 days after the onset of symptoms). 

Virus Detection and Isolation 

Viral nucleic acid was extracted from the swabs with the 
QIAamp Mini Viral RNA Extraction Kit and the DNeasy Blood 

6 Tissue Kit (Qiagen, Hilden, Germany). RT-PCR was 
performed with primers specific for CA24v and EV70 (CA24v- 
S: 5'-GTGAGTGCTTGCCCAGATTT-37CA24v-A: 5'- 
CTCCACTAGTGAGCGGTGTG-3'; and EV70-S: 5'-AGG- 
GATTCACCAGACATTGG-37EV70-A: 5'-ATTTTCCAC- 
CAGGCACTCTG-3') and nested PGR was used for the generic 
detection of adenoviruses [20]. The PGR products were analyzed 
with l"/o agarose electrophoresis. All ,52 swabs were cultured in 
fresh monolayers of HEP-2 cells. The cultures were incubated at 
36°C and observed daily for a cytopathic effect. Two blind 
passages were performed when no cytopathic effect was observed. 

Nucleotide Sequencing and Phylogenetic Analyses 

The 3Cpro and VPl genes of the isolates were amplified from 
the RNAs extracted from the culture supernatants and cycle 
sequenced with primers 3C1/3C2 and CA24v-2407S/CA24v- 
3438A, respectively [6] . The products were analyzed with an ABI 
PRISM 3100 Genetic Analyzer (Applied Biosystems, Hitachi, 
Japan). The MEGA 4.0 software [21] was used for the 
phylogenetic analysis. The homologies among the genotypes were 
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Nucleotide Sequence Accession Numbers 

All sequences reported in this study were deposited in GenBank 
under accession number JN788169-JN788308. In total, 276 
strains with complete 3Cpro sequences and 139 strains with 
complete VPl sequences that were available in GenBank (isolated 
between 1970 and 2011) were investigated. Sixty-three 3Cpro 
sequences and 41 VPl sequences, selected for their representa- 
tiveness, were used to construct the phylogenetic trees. Of these 
strains, 38 had complete sequen[:es for both the 3Cpro and Vpl 
regions. The GenBank accession numbers of the downloaded 
strains are shown on the trees (Fig. 2 and Fig. 3). 
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